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Methane- oxidizing bacteria (MOB) gain energy from the oxidation of methane and may play important 
roles in freshwater ecosystems. In this study, the community structure of planktonic MOB was investigated 
in a subtropical reservoir. Bacterial community structure was investigated through the analysis of the 16S 
rRNA gene. Three groups of phylogenetically distinct MOB were detected in the clone libraries of 
polymerase chain reaction products obtained with universal primers. The groups belonged to the class 
Gammaproteobacteria, the class Alphaproteobacteria, and the candidate phylum NCIO. The last group, 
which consists of close relatives of the nitrite reducer 'Candidatus Methylomirabilis oxyfera', was frequently 
detected in the clone libraries of deep-water environments. The presence of 3 groups of MOB in deep water 
was also shown by a cloning analysis of the pmoA gene encoding particulate methane monooxygenase. The 
dominance of 'M. oxyfera' -like organisms in deep water was confirmed by catalyzed reporter deposition- 
fluorescence in situ hybridization, in which cells stained with a specific probe accounted for 16% of total 
microbial cells. This is the first study to demonstrate that close relatives of the nitrite reducer can be major 
component of planktonic MOB community which may affect carbon flow in aquatic ecosystems. 

Methane- oxidizing bacteria (MOB) are capable of gaining energy from the oxidation of the simplest 
hydrocarbon, methane. This ability is observed only in restricted lineages of prokaryotes, and 4 major 
phylogenetic groups of MOB are currently known. Two of these MOB groups, belonging to the classes 
Gammaproteobacteria and Alphaproteobacteria in the phylum Proteobacteria, have been well characterized. 
These 2 groups are also referred to as type I and type II, respectively. Most type I MOB are members of the 
feLYnily Methylococcaceae, but the hmily Methylothermaceae, was recently proposed\ Type II MOB species belong 
to the families Methylocystaceae and Beijerinckiaceae. The third group, MOB in the phylum Verrucomicrobia, is 
extremely acidophilic, and the family Methylacidiphilaceae was proposed to encompass this groups. The fourth 
group has no representative isolated in pure culture. Members of this group are recognized as 'Candidatus 
Methylomirabilis oxyfera' and its close relatives, belonging to the candidate phylum NCIO^'^. 'M. oxyfera has 
the notable ability to produce oxygen from nitrite, which enables nitrite- dependent methane oxidation without an 
external supply of oxygen^. Accordingly, the dominance of 'M. oxyfera and its relatives has been observed 
consistently in methane -oxidizing enrichment cultures established under anoxic nitrite-reducing conditions'"^. 
Based on these studies, 'M. oxyfera' -like bacteria detected in natural environments have been regarded as the 
MOB responsible for nitrite- dependent methane oxidation^"^^. 

Methane oxidation by all known MOB is mediated by monooxygenases, in contrast to the reverse methano- 
genesis pathway used by methane- oxidizing archaea. Especially, the particulate methane monooxygenase is 
found in most of the known MOB. The pmoA gene encoding the a-subunit of this enzyme has been widely used 
as a marker to detect and identify MOB, along with the 16S rRNA gene. By using these and other markers such as 
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phospholipid fatty acids, the community structures of MOB have 
been intensively investigated in various types of ecosystems, mainly 
targeting type I and type II MOB. 

Considering their unique function, MOB may play important 
roles in the water columns of freshwater lake ecosystems. In the 
sediments of freshwater lakes, methanogenesis is a major terminal 
step in organic matter degradation. Methane release from sediments 
is consumed and assimilated by planktonic MOB, and therefore, they 
reduce methane emission from lakes and retrieve carbon released 
from the degradation of organic matter. The importance of MOB in 
the diets of other planktonic organisms has been discussed and eval- 
uated in previous studies^^"^^. The community structures of plank- 
tonic MOB have been investigated mainly in temperate ^^"^^ and 
subarctic^^'^^"^^ areas. These studies consistently indicated the dom- 
inance of type I MOB over type II MOB. Some studies detected type I 
MOB as the major components of entire bacterial communities in 
temperate water columns^ On the other hand, dominance of 
type II MOB in water column was reported in a tropical dam res- 
ervoir^^'^^. The primary factor to determine community structure has 
not been identified, although the involvement of temperature is 
likely. One possible approach to resolving this issue is to analyze 
the community structure of planktonic MOB in other climate zones 
where dominant MOB have not been identified. Especially, studies in 
subtropical zone may be valuable to fill a gap between temperate and 
tropical zones. 

In this study, the community structure of MOB inhabiting 
the water column of a subtropical reservoir was investigated to 
complement the current knowledge of planktonic MOB in fresh- 
water environments. The study site, the Feitsui Reservoir (FTR), is 
located in northern Taiwan^^'^°. The reservoir is protected from 
human activities to maintain water quality. 

Results 

Physicochemical factors. Vertical profiles of physicochemical 
factors measured in 4 sampling days are shown in Fig. 1. Water 
temperature ranged from 16.6 to 31.1°C. Depth-related changes in 
the availability of methane and oxygen were consistently observed. 
Oxygen was almost depleted at 90 m (2% or less) in all sampling 
days. A high concentration of methane was observed under anoxic 
conditions, and only a trace amount of methane was detected in the 
oxic water column. In Nov 2012, methane concentration was not 
determined on the day when samples for DNA-based analysis were 
obtained (Nov 13th), but measured on Oct 23rd and Nov 27th. On 
those days, methane concentration at 90 m was 2.9 and 7.1 |iM, 
respectively, and did not exceed 0.1 jiM at the other depths. 

Detection of MOB in bacterial 16S rRNA gene clone libraries. The 

phylogenetic affiliations of 344 clones from 7 libraries of surface- and 
deep-water samples were analyzed (phylum-level composition of the 
libraries is shown in Supplementary Table SI), and 68 clones were 
identified to be MOB as described below. As type I and type II MOB, 
18 clones of the family Methylococcaceae and 14 clones of the family 
Methylocystaceae were obtained, respectively (Fig. 2). In addition, 36 
clones from the deep-water samples were also identified as MOB, 
since they were closely related to 'M. oxyfera\ with sequence iden- 
tities of 96% (show as "NCIO" in Fig. 2). 

These MOB clones were further grouped into species-level OTUs 
(Supplementary Table S2). All clones of type II belonged to a single 
OTU (OTU_Al), but clones of type I were grouped into 8 OTUs 
(Fig. 3). The most frequently detected type I OTU was OTU_Gl, 
which included 10 clones (Supplementary Table S2). The 
'Methylomirabilis -like clones were very closely related to each other, 
and grouped into a single OTU, OTU_Nl. In fact, these clones could 
not be differentiated even at a cutoff value of 0.01. The OTU_Nl was 
the most frequently detected species-level phylotype among MOB in 
all deep-water libraries (Supplementary Table S2). An environmental 



clone from the sediment of a freshwater lake, clone B84 from Lake 
Biwa^\ had very high sequence identities with OTU_Nl clones 
(Fig. 4). The clone B84 was classified in the same OTU with these 
clones, in case that it was included in the grouping with the same 
methods. 

Cloning analysis of pmoA genes. From the 2 libraries constructed 
with the primer pair A189/mb661r, 101 clones were sequenced. They 
were grouped into 23 OTUs, 6 of which were detected in both 
libraries (Supplementary Fig. SI). Phylogenetic analysis revealed 
that only OTU_PP22 belonged to the type II lineage and that the 
other OTUs were relatives of type I MOB (Supplementary Fig. S2). 

From the libraries constructed with the primer pair A189/682_ 
NCIO, 11 clones were sequenced (4 from Nov 2012 at 90 m, 7 from 
Dec 2012 at 90 m). All clones sequenced were closely related to each 
other and to 'M. oxyfera. When the protein sequences of these clones 
were compared with those of 'M. oxyfera, only 5-7 amino acid 
differences were observed in 166 amino acid positions. Further, these 
clones were very closely related to 21 clones from the sediment of 
Lake Biwa, obtained with the same pair^\ In pairwise comparisons of 
these 32 clones (21 from FTR and 12 from L. Biwa), only 4 or fewer 
amino acid sequence differences were found. Phylogenetic position 
of a representative clone is shown in Fig. 5. 

Catalyzed reporter deposition-fluorescence in situ hybridization 
(CARD -FISH). The dominance of 'M. oxyfera -like: organisms in 
deep water suggested by DNA-based analyses was confirmed by 
CARD-FISH analysis. The water sample obtained at 90 m in Dec 
2013 contained cells stained with the specific probe DBACT-1027 at 
a concentration of 1.1 X 10^ cells ml"^ (Supplementary Fig. S3). The 
cells positive for DBACT-1027 comprised 16% of cells stained with 
DAPI and 21% of cells stained with EUB I-III, and were significantly 
frequently observed than cells hybridized with the probes Mg669 and 
Ma464 (p < 0.01, for both), targeting type I and type II MOB 
respectively (Fig. 6). 

Discussion 

In this study, the most predominant planktonic methane oxidizer 
was 'M. oxyfera -like bacterium belonging to NCIO in deep-water 
samples (Fig. 2, Fig. 6). Considering that 'M. oxyfera -like pmo A gene 
sequences were also detected from the same samples, it is reasonable 
to consider OTU_Nl as a methane oxidizer. At present, a mech- 
anism for the dominance of OTU_Nl in the FTR is almost comple- 
tely unknown. To our knowledge, such a high relative abundance of 
the 'M. oxyfera -like bacterium has never been reported from other 
aquatic environments. Further research is needed to determine 
whether the dominance of this lineage is specific to the study site 
or if it is also observed in other sites in similar environments. Water 
columns have not previously been considered a major habitat of 'M. 
oxyfera -like bacteria, but subtropical freshwater lakes may become 
important for studying the diversity and distribution of these organ- 
isms as planktonic MOB. 

The clone library analyses of 2 genes suggested that type I diversity 
in the FTR was greater than type II diversity. In these analyses, only 
1 OTU of type II was detected in the libraries of the 2 genes. 
Considering their phylogenetic positions, OTU_PP22 may corre- 
spond to OTU_Al belonging to the genus Methylocystis. For type 
I, the correspondence between the OTUs of the 16S rRNA and pmo A 
genes is not easy to identify. The most frequently detected type I 
OTU of the 16S rRNA gene was OTU_Gl, suggested to be a 
Methylomonas species (Fig. 3). The most frequently detected OTU 
of the pmo A gene, OTU_PP19, was also a close relative of a 
Methylomonas species (Supplementary Fig. SI, Fig. S2). These results 
suggest that Methylomonas species are major type I MOB in the FTR. 
In previous studies using methods similar to those used in this study. 
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Figure 1 | Vertical profiles of temperature, dissolved oxygen, methane, nitrate, and sulfate in FTR. 



Nov, 2012 



Dec, 2012 



10 m 



90 m 




Jul, 2013 

!■ Type I I 
'■ Type II I 
■ NC10 

IH other I 



Dec, 2013 



^ ^ » ^ 

67 n = 53^^ n = 48 n = 44 



Figure 2 | MOB in clone libraries of bacterial 16S rRNA gene. Total 
number of clones analyzed in each library is shown by "n". 



the dominance of Methylobacter species among planktonic MOB was 
shown in the water columns of freshwater environments 

In the CARD-FISH analysis of the water sample obtained at 90 m 
in Dec 2013, there was no significant difference between numbers of 
cells hybridized with Ma464 and Mg669 (Fig. 6). The probe Ma464, 
used to estimate type II MOB abundance, perfectly matches all type II 
clones obtained in this study. In contrast, the probe Mg669 had 1-2 
mismatches against 6 type I clones belonging to 5 OTUs. In the water 
sample analyzed with CARD-FISH, OTU_Gl and OTU_G2 were 
detected as a type I MOB in the clone library (Table S2). Although 
these OTUs matches the probe perfectly, the number of detected 
clones is obviously too small to exclude the presence of other OTUs 
with mismatches. Generally, the relative abundance of type I and 
type II was largely different among the libraries (Fig. 2), but their 
clone numbers in the respective libraries are too small to discuss 
quantitatively. Further experiments including quantification by 
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Figure 3 | Phylogenetic relationships of the OTUs of 16S rRNA gene clones belonging to the families Methylocystaceae and Methylococcaceae. 

Numbers on nodes are percentage values of 1000 bootstrap resampling (values larger than 50 are shown). 



CARD -FISH with other probes will complement the results of the 
present study. 

In freshwater lake sediments, methane is produced in the main 
terminal process of organic matter degradation. Therefore, plank- 
tonic MOB have been considered a key link between sediment and 
pelagic carbon flow in freshwater lakes. Their contributions to 
carbon dynamics have been assessed by isotopic signature, assum- 
ing that MOB cell material is composed of carbon originated from 
methane^^"^^. Considering some recent findings, however, this 
basal assumption may not always be applicable. Vermcomicrobial 
MOB are known to oxidize methane but do not assimilate it as a 
carbon source^^'^^. They fix inorganic carbon via the Calvin-Benson- 
Bassham cycle at the expense of energy derived from methane 



oxidation. In addition, recent research suggests that methane- 
dependent autotrophy is also a physiological feature of 'M. oxyfera'^'^. 
In the ecosystems dominated by these recently recognized MOB, the 
re-entry of methane to the organic carbon flow may be difficult to 
estimate. On the other hand, another carbon flow path via type I 
MOB was also suggested recently^^ In that report, a methane-oxid- 
izing bacterium was shown to excrete organic acids and hydrogen 
under oxygen-limited conditions, and genes responsible for this pro- 
cess are shared by various type I MOB^^. In the present study, type I 
MOB were detected in samples characterized with low oxygen con- 
centrations (Fig. 1, Fig. 2) and they were also detected in oxygen- 
depleted zones of other freshwater water columns^^'^^'^^. These MOB 
might play a particular role in carbon flow without the assimilation 



100 



0.005 



70| Uncultured bacterium clone 12-48 (AF351217) 



100 



j unc 



99 



100 l_f 
58 



Uncultured bacterium clone 12-24 (AF351214) 
"-Uncultured bacterium clone JMYB12-12 (FJ810544) 
Uncultured bacterium clone B84 (AB661586) 
0TU_N1 

Uncultured bacterium clone 49S1_2B_62 (DQ837259) 

Candidatus 'Methylomirabilis oxyfera' (NR_1 02979) 

Uncultured bacterium clone 49S1_2B_43 (DQ837250) 



-Uncultured bacterium clone UA_12 (JX1 20382) 

Uncultured bacterium clone BacC-s_046 (EU335192) 



Figure 4 | Phylogenetic position of a clone belonging to the OTU_Nl. Numbers on nodes are percentage values of 1000 bootstrap resampling (values 
larger than 50 are shown). 
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Figure 5 | Phylogenetic position of PmoA sequence deduced from a pmoA gene clone obtained in this study (12Nov_NP27), with the primer pair 
A189/682_NC10. Numbers on nodes are percentage values of 1000 bootstrap resampling (values larger than 50 are shown). 



of methane by converting it to carbon dioxide and organic acids 
available for other functional microbe groups. 

Methods 

Sample collection and chemical analyses. Samples were obtained from the Feitsui 
Reservoir (FTR) located in northern Taiwan (121° 34' E, 24° 54' N)^'. Water 
sampling was performed at the same site as that of a previous study^" using a 5-L Go- 
Flo bottle. In the present study, samples collected on the following 4 sampling days 
were analyzed: 13 Nov 2012, 11 Dec 2012, 30 Jul 2013, and 24 Dec 2013. Water 
samples were collected from 6 depths (0, 10, 30, 50, 70, and 90 m from the water 
surface) and were subjected to chemical analyses. Additional samples were obtained 
from a depth of 100 m in Dec 2013 and from a depth of 95 m in Jul 2013. 

Dissolved methane concentration was determined with headspace equilibration 
methods^*^'^^. Briefly, portions of the water samples were poured into glass vials. The 
vials were immediately sealed with butyl rubber stoppers leaving no headspace. The 
vials were stored in a cooling box (around 4°C) until addition of HgCl2 solution 
conducted within 3-4 h after the sampling. He gas (>99.999% purity) was injected 
into the bottles containing the fixed water to create headspace, and the vials were 
shaken vigorously to equilibrate. Methane in the headspace was quantified using a gas 
chromatograph (GC-2014; Shimadzu) equipped with a flame ionization detector to 
calculate dissolved methane concentration in the water. The concentrations of nitrate 
and sulfate were measured with an ion chromatograph (DX-120; Dionex), and dis- 
solved oxygen concentration was determined by Winkler titration. Exceptionally, in 
the sampling day of Dec 2012, dissolved oxygen was measured with an oxygen sensor 
calibrated by the titration. 

Polymerase chain reaction (PGR) of 16S rRNA genes. Planktonic microorganisms 
in the water samples were collected on 0.22- |am filters (Sterivex filter cartridges; 
Millipore). The filters were immersed in RNA-later (Ambion) and kept at 4°G until 
DNA extraction. DNA samples were prepared using the bead-beating method. Prior 
to DNA extraction, the filters were washed with Tris-EDTA (TE) buffer to remove 
RNA-later. The washed filters were cut into small pieces and then placed into screw- 
cap tubes containing 600 |il of TE buffer, 30 |il of 20% sodium dodecyl sulfate, and 
0.5 g of glass beads (equal- weight mixture of 0.1- and 0.05-mm-diameter beads). 
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Figure 6 | Abundance of cells hybridized with probes targeting bacteria, 
MOB of typel, typell, and NCIO. Error bars represent ± SD. 



After adding 600 [il of phenol-chloroform-isoamyl alcohol (25 : 24 : 1, v/v/v), the 
samples were shaken twice with a bead beater (FastPrep-24; MP Biomedicals) at 
4.0 m s"^ for 30 s, followed by centrifugal separation. DNA in the aqueous 
supernatant was precipitated by adding NaGl and isopropanol, and excess salt was 
removed by rinsing the pellet with 70% ethanol. 

From 7 DNA samples at 2 water depths (10 and 90 m on 4 sampling days, except 
10 m in Jul 2013), fragments of bacterial 16S rRNA genes were amplified using the 
universal primer pair 27F and 1492R. PGR amplification was initiated with 2 min of 
denaturation at 94°G, followed by 28 cycles of 30 s at 94°G, 30 s at 55°G, and 90 s at 
72°G. Final extension was carried out for 10 min at 72°G. Glone libraries were con- 
structed from the resulting PGR products from the 7 samples using the pGR2.1- 
TOPO vector and competent TOP 10 cells (Invitrogen). For clones randomly selected 
from the libraries, inserts of the 16S rRNA gene were fully sequenced. After exclusion 
of chimeras and other anomalous sequences using Mallard software^*^, all sequences 
were classified using the Ribosomal Database Project classifier to sort out MOB 
clones. MOB clones were grouped into operational taxonomic units (OTUs) using 
mothur software^'^ at a cutoff value of 0.03. The clone sequences representing OTUs 
were aligned with related sequences using the program GlustalX'*". Phylogenetic trees 
were constructed using the neighbor- joining method in the program MEGA5^\ 

Cloning analysis oipmoA genes. Fragments of the pmoA genes were amplified with 2 
primer pairs from 2 DNA samples of water obtained from a depth of 90 m in 2012. 
One of the primer pairs, A 189^^ and mb661r'*^ was used to detect proteobacterial 
MOB. PGR was initiated with 5 min of denaturation at 94°G, followed by 34 cycles of 
30 s at 94°G, 30 s at 56°G, and 45 s at 72°G. Final extension at 72°G was carried out 
for 10 min. The other primer pair, A189 and 682_NG10^\ was used to detect pmoA 
genes of 'Methylomirabilis -like organisms. PGR conditions for this primer pair were 
almost the same as those for A189/mb661r, but initial denaturation was shortened to 
2 min, and annealing temperature was 55° G. Four clone libraries (2 samples, 2 primer 
pairs) were constructed as described above. The resulting sequences were translated 
into amino acid sequences, and distance matrices were calculated with Poisson model 
using MEGA5. The distance matrices generated were imported into mothur to group 
the clones into OTUs at a cutoff value of 0.03. Phylogenetic trees were constructed 
with amino acid sequences of the clones representing respective OTU, using the 
neighbor-joining method. 

CARD-FISH. GARD-FISH analysis was performed on the water sample obtained 
from a depth of 90 m in Dec 2013. The lake water sample of lake water was mixed with 
4% paraformaldehyde solution in phosphate-buffered saline to obtain a final 
concentration of 0.9%. After fixation at 4° G for ca. 4 hour, fixed cells were trapped 
onto track-etched polycarbonate membrane filters (Gyclopore, 0.22-|am pore size, 25- 
mm diameter; Whatman) and stored at — 20°G until further processing. GARD-FISH 
analysis methods were similar to those described previously'^'^'^^ Prior to successive 
treatments in liquids, filters were treated with 0.1% (w/v) low- melting-point agarose 
(SeaKem, FMG BioProducts) to hold the cells. The filters were cut into small pieces, 
and one of them was used for direct count of total cell number by staining with 4',6- 
diamidino-2-phenylindole (DAPI). Other pieces for GARD-FISH analysis were 
incubated in lysozyme solution (10 |ig/mL lysozym in 0.1 M Tris-HGl, 0.05 M 
EDTA [pH 8.0]) for 25 min at 37°G to permeabilize the cells. To eliminate 
endogenous peroxidase activity, the filter pieces were treated with 0.01 M HGl with 
0.3% H2O2 for 10 min at room temperature. The samples were hybridized using 
probes labeled with horseradish peroxidase. The specific probe DBAGT-1027^ was 
used to detect 'Methylomirabilis -like organisms. It perfectly matches the sequences of 
NGIO from the FTR and does not match the other sequences obtained in this study. 
Probe mix EUB 1-llP^'^^ was used to estimate bacterial abundance. Two probes, 
Mg669 and Ma464'^^ were also used to detect type 1 and 11 MOB, respectively. 
Hybridization was performed for 2 h at 46°G. Formamide concentrations in the 
hybridization buffers were 20% (for Mg669 and Ma464), 35% (for EUB l-IIl), or 40% 
(for DBAGT-1027). After washing at 48°G, signal amplification was performed with 
Alexa Fluor 488 tyramide at 46°G for 15 min. GARD-FISH counts were calculated as 
means from 15-31 randomly chosen microscopic fields (on average, approximately 
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1200 stained cells were counted for each probe). Results from the group -specific 
probes were compared by Tukey- Kramer test. 

Nucleotide sequence accession numbers. Nucleotide sequences obtained in this 
study have been assigned the DDBJ/EMBL/GenBank accession numbers AB930495- 
AB930950. 
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